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We employ the chiral nucleon-nucleon potential derived using the method of unitary 



O ! transformation up to next-to-next-to-leading order (NNLO) to study bound and scatter- 



ing states in the two-nucleon system. The predicted partial wave phase shifts and mixing 
parameters for higher energies and higher angular momenta beyond the ones which are 
fitted are mostly well described for energies below 300 MeV. Various deuteron properties 
are discussed. We also apply the next-to-leading order (NLO) potential to 3N and 4N 
^ ! systems. The resulting 3N and 4N binding energies are in the same range what is found 
*^ ' using standard NN potentials. Experimental low-energy 3N scattering observables are 
also very well reproduced like for standard NN forces. Surprisingly the long standing 
^ ' ^y-puzzle is resolved at NLO. The cut-off dependence of the scattering observables is 
j^ ■ rather mild. 

1. INTRODUCTION 

Over the past decade effective field theory methods have been successfully applied to 
many different physical processes and became, in particular, a standard tool in calculating 
low-energy properties of the purely pionic and pion-nucleon systems. In this method one 
starts from the most general Lagrangian for pions and nucleons. Apart from the usual 
symmetry constraints, spontaneously and explicitly broken chiral symmetry of QCD is 
taken into account. As a consequence, the couplings of pion fields, which play the role 
of Nambu-Goldstone bosons within this formalism, are of derivative type. As has been 
shown by Weinberg more than 20 years ago, one can use such an effective Lagrangian to 
derive the expansion of low energy S-matrix elements in powers of some scale Q related 
to low external momenta of pions and nucleons JIl (the 4-momenta of the external pions 
and the 3-momenta of the external nucleons). 

Motivated by successful applications of CHPT in the vrvr and vrA^ sectors, Weinberg pro- 
posed in 1990 to extend the formalism to the NN interaction ^. The crucial difference 
is, however, that the nucleon-nucleon interaction is non-perturbative at low energies: it 
is strong enough to bind two nucleons in the deuteron. Thus, direct application of CHPT 



to the NN amplitude will necessarily fail. The way out of this problem, proposed by 
Weinberg, is to apply CHPT not to the amplitude but to a kernel of the correspond- 
ing integral Lippmann-Schwinger (LS) equation using time-ordered perturbation theory. 
Such a kernel or effective NN potential is then defined as a sum of all diagrams with two 
incoming and outgoing nucleons that do not contain purely nucleonic intermediate states 
(reducible diagrams). A systematic power counting for the potential can be derived in a 
similar manner as for the vrvr and irN systems. These ideas have been extensively studied 
by the Texas-Seattle group [Q. They have obtained an energy dependent two-nucleon 
potential^ at NNLO. With altogether 26 free parameters^] they were able to achieve a 
qualitative agreement with the NN data. 

Alternatively, one can also do the expansion directly on the level of the scattering am- 
plitude, as it has been proposed by Kaplan et al. 0. In their power counting scheme the 
leading order non-perturbative contribution in the S-channels is given by resummation of 
the iterated non-derivative NN contact interaction. All corrections to this result including 
those from various pion exchanges are treated as small perturbations 0. 

In what follows we will concentrate on the potential approach to the few-nucleon scat- 
tering. In ref. [^ (referred to as I from here on) we constructed the NN and 3N potential 
based on the most general chiral effective pion-nucleon Lagrangian using the method 
of unitary transformations. For that, we developed a power counting scheme consistent 
with this projection formalism. In contrast to previous results obtained in old-fashioned 
time-ordered perturbation theory, the method employed leads to energy-independent 
potentials. To LO, the potential consists of the undisputed OPE and two NN contact in- 
teractions. Corrections at NLO stem from the leading chiral TPE and additional contact 
terms with two derivatives. At NNLO one has to include the subleading TPE which con- 
tains pion-nucleon interactions with two derivatives (pion mass insertions). In contrast 
to ref. [0, we take the novel LECs from systematic studies of pion-nucleon scattering 
in CHPT [0]. The resulting effective potential is renormalized and applied to two-body 
bound state and scattering problems |^ (referred to as II). 

We also consider 3N and 4N systems within the chiral effective field theory |^ (referred 
to as III). To the best of our knowledge this is the first time that xPT has been practically 
applied to nuclear systems beyond A = 2 within the Hamiltonian approach. In this first 
application we restrict ourselves to the NLO NN potential. At this order no three-nucleon 
forces (3NFs) occur[| ||^. The NLO results presented here are therefore parameter free 
and can serve as a good testing ground for the usefulness of the approach. Of course, 
some aspects of the 3N system have already been studied in nuclear EFT |10,|ll|, but not 
as direct extensions of the NN system as done here. 

Our manuscript is organized as follows. In Section 2 we will introduce the effective NN 
potential up to NNLO and give explicit expressions in momentum space. In Section 3 we 
apply it to the NN bound state and scattering problems. In the next section 4 we discuss 
NLO results for 3N and 4N systems. Our findings are summarized in Section 5. 



^Such explicit energy dependence of the potential complicates its application to three- and more-nucleon 

problems. 

^Some of them are redundant and can be eliminated from the Lagrangian via Fierz transformation. 

•^The 3NF appears first at NNLO and contains 5 undetermined parameters. 



2. THE NN POTENTIAL 

In I, we derived the two-nucleon potential at next-to-leading order. To leading order, 
it consists of the one-pion-exchange potential (OPEP) and two S-wave (non-derivative) 
contact interactions. The latter are parametrized by two LECs Cs and Ct 0, 
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where /^ = 92.4 MeV is the pion decay constant, M^^ = 138.03 MeV the pion mass and 
Qa = 1-26 the axial-vector coupling. Further, q = p ' — p, where p and p ' stand for 
initial and final nucleon momenta, respectively. Clearly, the LO potential can only serve 
as a very crude approximation of the NN force. At NLO order, we have three distinct 
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Figure 1. NNLO corrections to the NN potential. The filled squares denote tttiNN 
vertices with two derivatives. 



contributions. The first one represents the one-loop self-energy and vertex corrections to 
the OPEP and to LO contact interactions. As discussed in II and |l^, such contributions 
do not lead to form factors and only renormalize the appropriate couplings. Secondly, 
one has seven different NN contact interactions with two derivatives: 
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AfAf,cont 



Ci g 2 + C2 P + (C3 g ' + C4 k^) (ai • 0^2) + tC^ l (^i + ^2) • (g x fc) 



+ Ce{q-ai){q-a2) + C7{k-a,){k-a2) . (2) 

Thirdly, there is the genuine TPEP. After performing the renormalization it reads: 
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L(g) = -v/4M2 + g2 1n ^m 



(4) 



Here P{k^ g ) represents the polynomial part. 



At NNLO we have the additional TPEP with one irnNN vertex with two derivatives as 
shown in Figure |1|. Three different interactions of this type enter the effective Lagrangian. 
For the corresponding LECs 01^3^4 we take the values obtained from fitting the invariant 
amplitudes inside the Mandelstam triangle, i.e. in the unphysical region [0. Various 
NNLO self-energy and vertex correction diagrams lead again only to renormalization of 
the corresponding couplings. Note further that the contact interactions are renormalized 
by an infinite amount in order to absorb all UV divergences of the NNLO TPEP. The 
explicit expressions for the renormalized NNLO TPEP can be found in II and in [|T^ . 

The effective potential discussed above is only meaningful for momenta below a certain 
scale and therefore needs regularization. That is done in the following way: 

Vip,p')^fAip)Vip,p')hip'), (5) 

where /a(p) is a regulator function. Here we use two different regulator functions: 

/f -P(p ) = ^(A^ - p') , fr'^^p) = exp(-pVA^) . (6) 

3. RESULTS FOR THE NN SYSTEM 

Having introduced the effective potential up to NNLO we will now apply it to the NN 
system. 

3.1. Phase shifts 

In order to apply the effective potential to bound state and scattering problems, one 
first needs to determine the unknown couplings related to the NN contact interactions. At 
LO one has two such constants (Cs, Ct), whereas nine constants (C5, Ct, Ci . . . C7) occur 
at NLO and NNLO, see eqs. (|l|), (H). To pin down this couplings we do not perform global 
fits as done in ref. 0. Rather we introduce independent new parameters by projections 
on to appropriate partial waves. To leading order, the two S-waves are depending on one 
parameter each. At NLO and NNLO, we have one additional parameter for ^Sq and ^5*1 
as well as one parameter in each of the four P-waves and in ei. We thus can fit each 
partial wave separately, which makes the fitting procedure not only extremely simple but 
also unique. All fits have been performed to the phase shifts of the Nijmegen phase shift 
analysis (NPSA) [0] for laboratory energies smaller than (50) 100 MeV at (NLO) NNLO 
as described in detail in II. Having determined the parameters, we can now predict the 
S- and P-waves for energies above 100 MeV. All other partial waves are parameter free 
predictions for all energies considered. In Figure |^ we show the results using the sharp 
regulator function with A = 500 and 875 MeV for LO (NLO) and NNLO, as discussed in 
II. At NNLO both S-waves agree well with the data. In the case of ^5*0 the OPEP plus 
non-derivative contact terms is insufficient to describe the phase shift (as it is well-known 
from effective range theory and previous studies in EFT approaches). This also holds true 
for some P-waves. It is known that TPEP alone is too strong in some of the D- and F- 
waves. To the order we are working, there are no contact interactions so that apart from 
varying the cut-off we have no freedom here. Still, for our best global fit, the D-waves 
are quite reasonably described. The result for the ^^3 phase is rather astonishing, since 
e.g. in the Bonn potential [|1^ the correlated two-pion exchange is very important in the 



description of D3. Such tttt correlations are only implicitly present at NNLO. They are 
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Figure 2. Predictions for various partial wave phase shifts (in degrees) for nucleon 
laboratory energies E\^\, below 300 MeV (0.3 GeV). The dotted, dashed and solid curves 
represent LO, NLO and NNLO results, in order. The squares depict the Nijmegen PSA 
results. 



hidden in the strengths of the LECs C\ and C3[]TH]. From the higher partial waves the ^Gs 
is the most interesting one. As in the case of ^-Ds, the OPEP alone is not sufficient to 
describe the phase shifts at energies above 100 MeV, as it could be expected for such a 
high value of angular momentum. The leading TPEP makes the situation even worse. 
Only inclusion of the subleading TPEP at NNLO allows for a very accurate and parameter 
freeQ description of the phase shift. 



Table 1 

Scattering lengths, effective range and shape parameters for the ^5*1 wave at NLO and 

NNLO compared to the NPSA. 



a [fni] 



r [fm] 



f 2 [fm^ 



f3 [fm^ 



f4 [fm'^] 



3^1 NLO 


5.434 


1.711 


0.075 


0.77 


-4.2 


35i NNLO 


5.424 


1.741 


0.046 


0.67 


-3.9 


^S^ NPSA 


5.420 


1.753 


0.040 


0.67 


-4.0 



We have also calculated the values of various effective range and shape parameters for 



*The cut-off dependence in this channel is rather weak. 



^Si. In Table |l] we compare our NLO and NNLO predictions for these quantities with the 
results from the Nijmegen PSA. As expected, in all cases the NNLO results go closer to 
the data than the NLO ones. Note further that in both cases one has the same number 
of free parameters. Therefore we consider the improvement as an indication of a good 
convergence of our expansion. 

3.2. Deuteron 

We now turn to the bound state properties. At NNLO (NLO), we consider an ex- 
ponential regulator with A = LOS (0.60) GeV, which reproduces the deuteron binding 
energy E^ = —2.224575(9) MeV within an accuracy of about one third of a permille 
(2.5 percent)^: E^ = -2.2238 (-2.1650) MeV. The asymptotic D/S ratio, called t], and 
the strength of the asymptotic wave function. As, as well as the root-mean-square mat- 
ter radius r^ are well described at NNLO (NLO): r] = 0.0245 (0.0248), As = 0.884 
(0.866) fm~^/^ and r^ = 1.967 (1.975) fni compared to experimentally measured values 
r) = 0.0256(4), As = 0.8846(16) fm^^/^ g^^d rj = 1.9671(6) fm, respectively. The D-state 
probability, which is not an observable, seems to be most sensitive to small variations 
in the cut-off. We observe a small deviation in the values for the quadrupole moment, 
Qd = 0.262 (0.266) fm^ compared to the experimental value Q^ = 0.2859(3) fm^, which 
is also typical for conventional potentials and possibly arises from the missing exchange 
current operators. 

4. RESULTS FOR THE 3N AND 4N SYSTEMS 

4.1. Bound states 

In this section we present the binding energies for 3N and 4N systems calculated from 
our NLO potential. As already pointed out above no 3NF appear at this order. The NLO 



Table 2 

Theoretical '^H and ^He binding energies for different cut-offs A compared to CD-Bonn 
predictions and to the experimental '^H binding energy and the Coulomb corrected ^He 
binding energy. 

NLO 540 NLO 560 NLO 580 NLO 600 CD Bonn Exp. 

Esn [MeV] -8.284 -8.091 -7.847 -7.564 -8.012 -8.48 

E4Hc [MeV] -28.03 -26.91 -25.55 -23.96 -27.05 -29.00 



results shown here are therefore parameter free predictions and can serve as a good testing 
ground for the usefulness of the approach. In order to investigate the cut-off dependence of 
3N and 4N observables we have generated several NN potentials corresponding to different 
exponential cut-offs between A =540 and 600 MeV. They were all fitted to the ^5*0, ^5*1- 
^Di, ^Pi and ^Po,i,2 NN phase shifts up to Eiab = 100 MeV, see III for more details. We 
find for the fully converged solutions of the corresponding Faddeev-Yakubovsky equations 



^Note that the deuteron binding energy is not used to fit the free parameters in the potential. 



the 3N and 4N binding energies as given in Table 0. The ranges are compatible with what 
is found using realistic potentials. Note that the NN forces were included up to the total 
NN angular momentum oi jmax=Q- 

4.2. nd scattering 

Also, 3N scattering can be solved rigorously nowadays |]16| and we show in Figs. |]some 
elastic observables as a particular example, namely the angular distributions, the nucleon 
analyzing power Ay and the tensor analyzing power T2o. Except for Ay there are no nd 
data for those energies. The discrepancies between data and theory for T2k and for the 



differential cross section can be traced back to the pp Coulomb force [19|. Thus except for 
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Figure 3. nd elastic scattering observables at En = 3MeV (upper row) and En = lOMeV 
(lower row) for the chiral forces (A = 540 MeV/c, dotted curve; A = 600 MeV/c, dashed 
curve), and CD-Bonn (thick solid curve), pd data, are (+) [|l^. nc? data are (o) [|18| . 



Ay the agreement of CD-Bonn with the data is rather good, which is a well known fact 
and is just given for the sake of orientation. Our parameter free predictions for da and 
T2k agree rather well with the CD-Bonn result and thus with the data. We consider this 
to be an important result, demonstrating that the chiral NN forces are very well suited 
to also describe low-energy 3N scattering observables rather quantitatively. On top of 
that, surprisingly for us, the chiral force predictions are now significantly higher in the 
maxima of Ay than for CD-Bonn and break the long standing situation, that all standard 
realistic NN forces up to now underpredict the maxima by about 30 %. That used to be 
called the 74y-puzzle|]l6|. We are now in fact rather close to the experimental nd values. 
Note that on a 2N level the chiral potential predictions for the analyzing power agree 
well with the predictions based on the Nijmegen PSA |^. Our predictions for various 
break up observables agree also very well with the ones obtained from CD-Bonn, as it is 



demonstrated in 



5. CONCLUSIONS 



Our findings show that the scheme originally proposed by Weinberg works quantita- 
tively and much better than it was expected. It extends the successfull applications of 
effective field theory (chiral perturbation theory) in the pion and pion-nucleon sectors to 
systems with more than one nucleon. Most of the np partial waves are well described at 
NNLO. The NNLO TPEP turnes out to be too strong in the triplet F-waves. This is 
expected to be cured at N^LO. 

The very first results using chiral NN forces in 3N and 4N systems are very promising 
and show that these effective chiral forces are very well suited to describe also low energy 
properties of nuclear systems beyond A = 2. They agree rather well with standard 
nuclear force predictions as exemplified with CD-Bonn and most importantly they break 
the stagnation in the A^^-puzzle. It will be very interesting to perform the next step and 
use the NNLO NN forces, which would also require inclusion of the leading 3NF. 
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